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We have derived the explicit velocity distribution function of the steady-state Boltzmann 
equation for hard-core molecules to second order in density and the temperature gradient. 
We correct some mistakes in the results in this addenda. 

The second-order coefficients Bj} r in eq.(38) in our recent paper should be written in 
the form: 
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Values for the constants b 2r and bg are summarized in Table HI a, n, T are the diameter of 
hard-core molecules, density and temperature of gases, respectively. Owing to the properties 
of aespWical—c tactions , B$ n can be obtained by replacing 2(d z T) 2 - (d x T) 2 - 
{dyT) 2 and 2d 2 T - d 2 x T - d 2 T in eq.© by Qd z Td x T and 6d z d x T, respectively, using an 
axis change. Similarly, C^J U can be obtained by replacing 2(d z T) 2 — (d x T) 2 — (d y T) 2 and 
2d 2 z T-d 2 x T-d 2 y T by 6d z Td y T and 6d z d y T, respectively. sg )n can be obtained by replacing 
2(d z T) 2 - (d x T) 2 - (d y T) 2 and 2d 2 z T - 3 2 X T - d 2 T by 6[(d x T) 2 - (d y T) 2 ] and 6(d 2 x T - 
d 2 T), respectively, via an axis change; by replacement by \2d x Td y T and 12d x d y T, 

respectively. 

The calculation of B^ r in Appendix D of our paper [lj should be also changed as follows. 
Equations (D-4), (D-5) and (D-6) become 
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and 
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respectively. Here m and k are the mass of hard-core molecules and the Boltzmann constant. 
No other corrections appear in the calculation of B^ r in Appendix D of our paper Q. We 
emphasize that the final form and the numerical values bi r , bn r and b 2r of the velocity 
distribution function to second order in eq.(39) of our paper |l[ which we have used to 
calculate thermodynamic quantities are unchanged. We also note that the 4th approximation 
620 in Table III of our paper [lj still deviates from Burnett's 4th approximation 620 by a factor 
1.003. |2| This deviation is considered to be due to errors in Burnett's calculation, as was 
indicated in our paper [lj]. Furthermore, we note a miss print in eq. (A-4) of our recent 
paper [lj]. Its correct form is 



k + 1 3 k 

Gk,r = 2fc + g {(fc + r + -)-Bfc+i,r - B k+1>r _ x } + 2fc _ ^ { B k-i,r ~ (r + 1)-B fc _ 1>r+1 }. 



(5) 



Finally, we correct the numerical values At s in each component of the kinetic temperature 
T[ for i = x, y and z. In eq.(44) of our paper 1], the values of should be identical with 
those of Ap in Table IV of the paper [jj. We also indicate that the second-order pressure 
tensor in eq.(43) of our paper |l( should be uniform from the solubility conditions for the 



third-order solution ri 



(3) 



[1] Kim H.-D. and H. Hayakawa: J. Phys. Soc. Jpn. 72 (2003) 1904. 
[2] D. Burnett: Proc. Lond. Math. Soc. 40 (1935) 382. 
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TABLE I: Numerical constants b^ r (upper) and 6^, (lower) in eq.(^Q) 
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Abstract 

The velocity distribution function of the steady-state Boltzmann equation for hard-core molecules 
in the presence of a temperature gradient has been obtained explicitly to second order in density 
and the temperature gradient. Some thermodynamical quantities are calculated from the velocity 
distribution function for hard-core molecules and compared with those for Maxwell molecules and 
the steady-state Bhatnagar-Gross-Krook(BGK) equation. We have found qualitative differences 
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confirmed that the steady-state BGK equation belongs to the same universality class as Maxwell 
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I. INTRODUCTION 



The kinetic theory has long history and various investigations have been carried out 



analytically numerically and experimentally. 



equation 
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Text books on the Boltzmann 



111] show the fact that the Boltzmann equation has attracted much 
interest among researchers of kinetic phenomena. It is well accepted that the Boltzmann 
equation is one of the most reliable kinetic models for describing nonequilibrium phenomena. 

A number of studies on the Boltzmann equation have been based on two representative 
,„ode l8 of m olecu,es: h ard -co r e m olecu,es and Maxwell m o.ecu,esE fl Maxwell 
molecules interact via a force that is inversely proportional to the fifth power of the distance 
r, that is, Maxwell molecules interact via a central potential proportional to r -4 . It is well 



known that calculations involving the Boltzmann equation 
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'or Maxwel 
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molecules become 
In order to obtain 



much easier than those for hard-core molecules. Q, 0, 
the velocity distribution function of a dilute gas system in a nonequilibrium state, various 

;hod and the Hilbert method, 



methods, such as the Chapman-Enskog method, the Grad met 
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have been presented. 
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201 ] They are well 



known as methods to derive the normal solutions of the Boltzmann equation. 

Other kinetic models besides the Boltzmann equation, such as the linearized Boltzmann 
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27 



28, 



and the Bhatnagar-Gross-Krook (BGK) equation 

n 

, |33J, have also been proposed in order to avoid mathematical difficulties 
in dealing with the collision term of the Boltzmann equation. It has been believed that 
results of those kinetic models approximately agree with those of the Boltzmann equation, 
though they are accepted as quantitatively different ones from the Boltzmann equation: 
e.g. the Prandtl number is 2/3 for the Boltzmann equation, while it is 1 for the BGK 



equation. 
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The BGK equation retains important properties of the Boltzmann equation, 
such as the conservation laws and the H-theorem. It has been felt that characteristics 
of molecules, such as hard-core molecules and Maxwell molecules, can be absorbed into 
the relaxation time r. The most important contribution concerning the solution of the 



BGK equation has been made by Santos and his coworkers. 
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Q, 28, 29,0, 31,Q,Q 



They have solved the steady-state BGK equation exactly and compared its exact solution 



with the Chapman-Enskog solution derived to arbitrary order. 
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28 



29 



30j They 



have concluded that the Chapman-Enskog solution is asymptotically correct, although the 



2 



Chapman-Enskog series diverge. 



26, 



22, Q, 29, 



They have also compared high order 



velocity distribution function for the steady-state BGK equation with that for the Boltzmann 
equation for Maxwell molecules under uniform shear flow. [31L I32I 33^ 
On the other hand, it had been 

^B»4& the — second 
order solution of the Boltzmann equation by the Chapman-Enskog method. [l| Importance 
of the second-order coefficients has been demonstrated for descriptions of shock wave profiles 



34 



35, 



36 



However, we have realized that Burnett's 



and sou nd p ropagation phenomena, 
solution 151 is not complete. Though there are various studies on the transport coefficients 
of the Bol tZm ann Ration to second orderflQQQQQQ o r even to third orderQ, 
we have found that nobody has derived the explicit velocity distribution function of the 
Boltzmann equation for hard-core molecules to second order. This is a result of mathe- 
matical difficulties, as was indicated by Fort and Cukrowski |40(] . For Maxwell molecules, 
Schamberg 3] has derived the precise velocity distribution function of the Boltzmann equa- 
tion to second order by the Chapman-Enskog method, while Shavaliev JJZj has derived it 
implicitly by the moment method. 

In this paper, we derive the explicit velocity distribution function of the steady-state 
Boltzmann equation for hard-core molecules to second order in density and temperature 
gradients by the Chapman-Enskog method. The derivation of it is of some physical signifi- 
cance. For example, it has been required for calculating a nonequilibrium effect on the rate 
of chemical reactionj^J, because the nonequilibrium effect on the rate of chemical reaction 
does not appear to first order and it cannot be calculated by the moment method. [4]] There 
also exists a need to confirm the existence of universal nonlinear nonequilibrium statistical 
mechanics from the microscopic viewpoint as mentioned in ref . , so that it is necessary to 
derive the explicit velocity distribution function of the steady-state Boltzmann equation for 
hard-core molecules to second order and compare it with that for Maxwell molecules. Addi- 
tionally, the derivation of the explicit second-order solution of the steady-state Boltzmann 
equation for hard-core molecules by the Chapman-Enskog method will contribute to an un- 
derstanding of the difference or the correspondence between the Chapman-Enskog method 
and the Grad method by the direct comparison of the solution by the Chapman-Enskog 
method with that by the Grad method. In the present paper, we also discuss the relation 
between the steady-state Boltzmann equation and the steady-state BGK equation which 
has not been fully understood yet. 
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The organization of this paper is as follows. In £|TTJ we will introduce the Chapman- 
Enskog method to solve the steady state Boltzmann equation. In £11111 we will derive the 
explicit form of the velocity distribution function of the steady-state Boltzmann equation for 
hard-core molecules to second order in the density and temperature gradients. In particular, 
the result for the first-order coefficients is shown in eq.(|33|). and the results for the second- 
order coefficients are shown in eqs.([37|) and ()38|). The velocity distribution function of the 
steady-state Boltzmann equation to second order for hard-core molecules is shown explicitly 
in eq. ()39|) and graphically in Fig|21 and compared directly with that for Maxwell molecules 
()4())1 in FigOl In §IV| we will apply the velocity distribution function to second order for hard- 
core molecules (|39|) and those for Maxwell molecules (|4Up and the steady-state BGK equation 
to a nonequilibrium steady-state system under steady heat conduction and calculate some 
thermodynamical quantities. We stress the existence of the qualitative differences among 
hard-core molecules, Maxwell molecules and the steady-state BGK equation in the pressure 
tensor (J4*3|) and the kinetic temperature (}4*4"|) . Our discussion and conclusion are written in 
fVlandlVH respectively. 

II. METHOD FOR SOLVING THE STEADY-STATE BOLTZMANN EQUATION 

Let us introduce the Chapman- Enskog method to solve the steady-state Boltzmann equa- 
tion in this section. Assume that we have a system of dilute gases in a steady state, with 
velocity distribution function f\ = f(r, vi). The appropriate steady-state Boltzmann equa- 
tion is 



with f[ = f(r,v[) and f 2 = f(r,v' 2 ): v[ and V 2 are postcollisional velocities of vi and v 2 , 
respectively. The relative velocity of two molecules before and after the interaction has the 
same magnitude g = |vi— V2I ; the angle between the directions of the relative velocity before 
and after the interaction is represented by %■ The relative position of the two molecules is 
represented by b, called the impact parameter, and e represents the orientation of the plane 
in which g and g' = — v' 2 belong. The impact parameter b depends on the kind of 



vi-V/i = J(/i,/ 2 ), 



(1) 



where the collision integral J(fi, f?) is expressed as 
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interactions between molecules. Note that x can De expressed as a function of b for a central 
force, (see FigJTJ 

Suppose that the velocity distribution function fi can be expanded as: 

h = fl 0) + A (1) + A (2) + • ■ ■ = A (0) (i + } + # + •■■). (3) 

fi°^ is the local Maxwellian distribution function, written as 



( m \ 2 



/l(0) = " W {2^) ' 6XP 



2 



(4) 



2nT(r 

with m mass of the molecules and k the Boltzmann constant. n(r) and T(r) will be identified 
later as the density and the temperature at position r, respectively. Substituting eq.(jHJ) into 
the steady-state Boltzmann equation (Q), we arrive at the following set of equations which 
we will solve completely in this paper: 

L[f[ W=^^fi°\ (5) 

to first order and 

^]0? ) = vi-V/«-J(/« I /«) I (6) 

to second order. The linear integral operator is defined as 

L[A (0) ]X! = J J f ff ) fi G \X' l -X 1 + X' 2 -X 2 )gbdMedv 2 . (7) 

The solubility conditions of the integral equation (J5J) are given by 

$ iVl • VjfWa = 0, (8) 

where $i is one of the collisional invariants: 

$1 = 1, $ 2 = mvi, $3 = ^rnw\. (9) 

Substituting eq.(jlj) into the solubility conditions (jBJ), it is seen that nnT is uniform in the 
steady state. We use this result in our calculation. Similarly, the solubility conditions of the 
integral equation © are given by 

$ iVl ■ VA (1) d Vl = 0. (10) 



To construct solutions of the integral equations (0) and © definite, five further conditions 
must be specified; we identify the density: 

n(r) = [ /idv! = / /fdv!, (11) 



the temperature: 

Me^>) s / n^Adv, = / ^i/r'dv, (12) 

and the mean flow: 

C = J mvi/idvi = y mvi/^dvi. (13) 

Here we assume that no mean flow, i.e. Co = 0, exists in the system. The introduction of 
these conditions distinguishes the Chapman-Enskog adopted here from the Hilbert method 
in which the conserved quantities are also expanded. P] We assert that conditions (fTTj) . (fT2^) 
and (fT3~|) do not affect all our results in this paper. It should be noted that, to solve the 
integral equations (0) and (jHJ), we should consider only the case in which the right-hand 
sides of eqs.(J3} and © are not zero: if the right-hand sides of eqs.© and © are zero, 
the integral equations (0) and © become homogeneous equations which do not have any 
particular solutions. [2] 

III. BURNETT'S METHOD 

A. A general form of the velocity distribution function 



To solve the integral equations (jHJ and (jEJ), Burnett |15j assumed a general form of the 
velocity distribution function: 

OO _ OO oo „ 

h = fi 0) E r!r > + 2 )BorS l (c?) + £ (2^) 7 £ r!r(A; + r + 2 )Ykr{Ci)sr ^ (c?)] • (14) 

r=0 " fc=l r=0 

Here Ci = (m/2fi:T) 1 / 2 v 1 is the scaled velocity and T(X) represents the Gamma function. 
S%(X) is a Sonine polynomial, defined by 

00 

(1 - ui)- k - 1 e-^ = ^r(p + H l)S p k (X)u p . (15) 

p=0 

Ykr{ c i) is a linear combination of spherical harmonic functions: 

Y kr (ci) = B kr Y k ( Cl ) + 2 £ |zgJ[BWyW( Cl ) + c£>Z«( Cl )], (16) 
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where Bk r , Bj^J and Cj£} are coefficients to be determined. Introducing the normal spherical 
coordinate representation for c 1; i.e. c\ x = c\ sin 6 cos 0, c\ y = c\ sin 6 sin (f> and c\ z = c\ cos 9, 
the spherical harmonic functions Y k (ci) : (ci) and Z^\ci) are expressed as 



2kT 



k 
2 



and 



and 



n(ci)= ( — ) c*P fe (cos0), (17) 



yW( Cl ) = (-If ( ^ ) ' ctPi p) (cos^) cosp0, (18) 



4 P) (ci) = ("l) p f — ) ' cfP^cos*) sin P 0, (19) 



m 

with the Legendre polynomial Pk(cos8) and the associated Legendre polynomial Pjf\cos6). 

Assumption of the velocity distribution function of eq.flHJ) has some mathematical ad- 
vantages in our calculation. Firstly, it is sufficient to determine the coefficients B^, because 
Bfa! and C^J can be always determined from B kr by a transformation of axes owing to the 
properties of the spherical harmonic functions (|T7|). (fTHj) and (fHJj) . We call B kr , B^J and C^J 
the family of B kr . Secondly, some important physical quantities are related to coefficients 
Bkr, Bfa! ,Cj£: e.g. the density (fTTjl . the temperature ((E)) and the mean flow (JTHJ) with fi 
in eq.((T3J) lead to the five equivalent conditions [lfl llfil] : 

-Boo = 1, Bio = B$ = clo* = 0, B x = 0. (20) 

Similarly, the pressure tensor P^ defined by 

/ 2k,T\z f°° 

Pij = [~^~J J dcimciiCij/i, (21) 

for i, i = , and , is related to the family of B 20 , „h,ch is the on ly family in wh,eh Bn r nett 
was interested. |JJ 

The coefficients B kr except for those in eq. (|2T)|) can be calculated as follows. Multiplying 
the steady-state Boltzmann equation (JTJ) by 

Qkr^i) = (k + i)v / ^(^)^ fe (c 1 )^ + ,(c?), (22) 



7 



and then integrating over (2KT/m) 1 / 2 c 1 , it is found that 



2kT\ 2 



m 



< Ci • VQkr > av +V 



2kT\ 2 



m 



< c lQkr > < 



[2kT\ 



III J(Qkr-Qkr)fif2gbdbdedc 2 d Cl , (23) 



\ m ) 

where < X > av represents [2kT /mf'/ 2 J X/idci, and Q' kr represents the postcollisional 
Qkr- We should calculate both sides of eq.(j2HJ) for every k and r, because eq.(J2HJ) leads to 
equations to determine Bk r , as will be shown in Appendices IU1 and ID1 

For convenience, we introduce Qkr and Ak r as the left-hand and right-hand sides of 
eq. (j2~3*j) . respectively, i.e. 



q 



kr 



2kT\ 2 

J < Cl • VQkr >av +V 



m 



2kT\ 5 



< CiQ kr > 

( 



(24) 



and 



A 



kr 



(^") /// J(Qkr-Qkr)fif2gbdbdedc 2 d Cl . (25) 



We will calculate Qkr and A^ r separately. The result of Qk r becomes 

Qkr 



_( )?[(r + -)(D h , r d x T + E k , r d y T + G k ,rd z T) 



Tm /Ly 2 
— {Dk,r-id x T + Ek )r -idyT + Gk,r-id z T)\ 



+ d x 



,2k7\ i 

n )*D k . 

m 



rn 



k:r 



+ d z 



n( ) 2 G fc ,r 

m 



(26) 



where diX denotes dX/di for i = x, y and z. Dk >r , Ek, r and Gk, r are functions of the family 
of Bk r , as is written in Appendix IA"1 



B. The collision term A& r 

Next, we calculate the collision term A^ r in eq. (|25|) . We should specify the kind of 
interactions of molecules so as to perform the calculation of the collision term Ak r ] the 
impact parameter b is explicitly determined by specifying the type of interaction. For hard- 
core molecules, the impact parameter b is given by the relation 

b = a cos ^, (27) 



S 



where a is the hard-core molecular diameter and x is the scattering angle (see FigJIJ. There- 
fore, Akr for hard-core molecules, i.e. A^ r , becomes jis| 

Ag. = y /VfcUx) - ^(0)] sin | cos |d X , (28) 
where F kr (x) is defined as 

F£(x) = III Q' fcr / 1 / 2 < ? 'Medc 2 dc 1) (29) 

and we have used F&(0) = (2«;T/m) 3 J J J Q kr f 1 f 2 g^dedc 2 dc 1 . Note that x = if b > a. 
For Maxwell molecules, the impact parameter b is given by the relation 

bdb = -H( X )d X , (30) 



12 



13 



14J Therefore, A kr for Maxwell molecules, i.e. Aj^, 



where #(x) is a function of X-|£> 
becomes 

4£ = r [F° kr ( X ) - F° kr (0)Mx)d X - (31) 

Since A^. has been calculated by Schambergyjj], we will calculate only Aj? r in this paper. 
Schamberg's result for A^. is briefly summarized in our web page. (2] 

From eq.pS)). it is sufficient to calculate F kr (x) f° r ^kr- The details of F kr (x) are written 
in Appendix [B] Several explicit forms of A^ r are also demonstrated in Appendices O and 
iDl From the definitions (|24j) and (|25|) . both sides of eq. ()23j) for arbitrary k and r can be 
calculated for hard-core molecules via 

«£■ = Ag., (32) 

which produces a set of simultaneous equations determining the coefficients B kr , as is ex- 
plained in Appendices O and |D| Here fi^ r denotes fl k r in eq. (j2T?j) for hard-core molecules. 



C. Determination of B kr 

We will determine the first-order coefficients B\ r and the second-order coefficients B\ l r 
in accordance with the previous two subsections, which corresponds to solving the integral 
equations © and (JEJ), respectively. Here the upper suffices I and II are introduced to specify 
the order of B kr . We have provided an example of our Mathematica program for calculating 



these coefficients. 



42] 



1. The first order 



We show the results of the first-order coefficients B\ r of which the solution of the integral 
equation (j5J), <f>i \ is composed. They can be written in the form: 

B - = s ^4vB^- (33) 



Values of the constants b\ r are given in Table |U The calculation of B\ r is explained in 
Appendix [0 It is seen that B\ r is of the order of the Knudsen number K, which means that 
the mean free path of molecules should be much less than the characteristic length for changes 
in the macroscopic variables. Though Burnett derived B\ r only to 4th approximation, i.e. 
B\ r for r < 4, we have obtained B kr for r < 7 in this paper. This ensures the convergence of 
all the physical quantities calculated in this paper. It should be mentioned that our values 
of B\ r for r < 4 agree with Burnett's values 15 1. Once B\ r have been calculated, -B^ 1 can 
be written down directly by replacing d z T by d x T by symmetry, owing to the properties of 
the spherical harmonic function Cj^ 1 can also be obtained similarly by replacing d z T 
by d y T. Note that other terms, such as BfJ 1 , do not appear because p in B^J 1 must be k 
or less from eq.(fTfij). Substituting all the first-order coefficients derived here into eq.lflfij). we 
can finally obtain the first-order velocity distribution function f^. 



(2) 

2. Solubility conditions for (f>\ 

Since the first-order velocity distribution function f± has been obtained, the solubility 
conditions of the integral equation © should be considered before we attempt to derive an 
expression for <ffi. The solubility conditions for <j)f\ that is, eqs.(|TU]) lead to the condition 

V • J (1) = 0, (34) 

where J^ 1 -*, i.e. the heat flux for , can be obtained as 



. m J J-oo 2 

= - 6 4(-)*4vT, (35) 
64 \TTm J o l 

with the appropriate value for bn listed in Table |l] It must be emphasized that, since 3^ 2 \ 

(2) 

i.e. the heat flux for /{ , does not appear, the solubility conditions of the steady-state 
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(2) 

Boltzmann equation for <j)\ lead to the heat flux being constant to second order. From 
eqs.(jnU) an d (G3), we also obtain an important relation between (VT) 2 and V 2 T, namely 

^£ + V 2 T = 0. (36) 
Owing to the relation (|36|). terms of V 2 T can be replaced by terms of (VT) 2 . 



3. The second order 

We write down the results of the second-order coefficients B l £ r of which <j>^ is composed. 
We can determine the second-order coefficients B^ r appearing in eq.fJHJ) as 

<=J^(VTf. (37, 

Values for the constants bo r are summarized in Table |HJ The calculation of B^ r is shown in 
Appendix[D] It must be emphasized that Burnett did not obtain the second-order coefficients 
B^ r because he was only interested in the family of the coefficient B^ which is related 
to the pressure tensor -Pjj.Q] The necessity of the second-order coefficients B^ r for the 
calculation of some physical quantities will be discussed later. We have calculated B^ r to 
7th approximation, i.e. B^ r for r < 6 in this paper. 

The other second-order coefficients B^ r in eq.(|16)l can be written in the final form: 



Bl = [2 {d z T? - {d x Tf - {d y Tf] . (38) 



Values for the constants &2r are summarized in Table ITTT1 The calculation of B^ r is explained 
in Appendix \D\ Our 4th approximation B^ Q deviates from Burnett's 4th approximation B^ 
by a factor 1.003. This deviation is considered to be due to errors in Burnett's calculation 
of the second and the third terms on the right-hand side of eq. (ID 16)1 . Although Burnett 
obtained B^ r only to 4th approximation, i.e. B^, for r < 3, we have obtained B l ^ r to 7th 
approximation, i.e. B^} r for r < 6 in the present paper. Owing to the properties of the 
spherical harmonic function 1^. B^! U can be obtained by replacing 2(d z T) 2 — (d x T) 2 — 
(d y T) 2 in eq.()38j) by Qd z Td x T using an axis change; similarly, C^J 11 can be obtained by 
replacing 2(d z T) 2 - (d x T) 2 - (d y T) 2 by 6d z Td y T, and Bj£ )n can be obtained by replacing 
2(d z T) 2 - (d x T) 2 - (d y T) 2 in eq.flSHD by 6[(d x T) 2 - (d y T) 2 } via an axis change; c£ )n by 
replacement by 12d x Td y T. Note that other terms, e.g. B^ 11 , do not appear because p in 
B^ 11 must be k or less, from eq.(fIH)l. 

11 



One can see that both of Bg. and Bf r are of the order of K 2 . As is mentioned in Appendix 
iDl we have found the fact that Bj} r for k = 4, 6 and 8 do not appear, which strongly suggests 
that Bj} r for all k greater than 2 do not appear. Thus, we expect Bj} r = for k ^ 2 although 
Burnett 15] had believed that they would appear. We have calculated the second-order 

(2) 

coefficients now, so that we finally obtain /{ by substituting the second-order coefficients 
obtained here into eqs.()14p and (fTBj) . 



4- The velocity distribution functions to second order 

The velocity distribution function for hard-core molecules which we have derived in 
flIII C II and IIII C 31 valid to second order is now applied to a nonequilibrium steady-state 
system under the temperature gradient along x-axis. It can be written as 

4096mJ 2 rV -^ 3. _,_ . 9 . 

+ 56266f^E r!6 » r(r+ 2 )S l (c) 

11 r>2 

+ ^r!6 2r (2c 2 -c2-ar(r + ^5(c 2 )]}, (39) 

r>0 

where the specific values for 6 lr , 6 r and b 2r are found in Tables IH ITT1 and 111 11 respectively, 
and J x corresponds to the x component of the heat flux in eq.()35|). Note that we have 
changed Ci to c. Figure 121 provides the explicit form of the second-order velocity distribution 
function of the steady-state Boltzmann equation <p^ 2 \c) for hard-core molecules with 7th 
approximation b 0r and b 2r , scaled by mJ 2 /n 2 K 3 T 3 . Figure El shows that the scaled <p^ 2 \c) 
for hard-core molecules is strained symmetrically. 

In order to compare the macroscopic quantities for hard-core molecules with those for 
Maxwell molecules, we also adopt the velocity distribution function for Maxwell molecules 
derived by Schamberg^^] for our calculation of the macroscopic quantities. We apply Scham- 
berg's velocity distribution function to a nonequilibrium steady-state system under the tem- 
perature gradient along x-axis, that is, we take out the differential coefficients with respect to 
time t from Schamberg's velocity distribution function and also use the solubility condition 
(|3~4*|) for Maxwell molecules. The precise velocity distribution function of the steady-state 
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Boltzmann equation for Maxwell molecules to second order finally becomes 

x xfoUi / m \ — t r/ Aoi/ 2\ , 4096mJ^ ^ , 3 2 

r=2,3 

+ £ r\b 2r (2cl -c\- c\)T(r + 7 -)Sl{c 2 )}}. (40) 

r=l,2 

As can be seen from eqs.(|39jl and (|40|). the explicit form of the velocity distribution function 
for hard-core molecules becomes the sum of an infinite series of Sonine polynomials, while 
the precise form of the velocity distribution function for Maxwell molecules is the sum 
of five Sonine polynomials. Figure El gives the direct comparison of the scaled <f>^s in 
eq. ()39|) with 4th, 5th, 6th and 7th approximation bo r s and &2rS with the scaled (p^ in 
eq.(Jinj). It should be mentioned that, as Fig|3] shows, the scaled (j)^ in eq.(j3HJ) has not 
yet converged to 4th approximation. Furthermore, in this paper, we adopt the precise 
expression of the corresponding velocity distribution function for the steady-state BGK 
equation, though we do not write down the explicit form of it, which can be reduced from 
the general form of the the Chapman- Enskog solution for the steady-state BGK equation to 



arbitrary order 26J, |27|, |2a, |29j, |30 ] . We mention that for the steady-state BGK equation 



is identical to that for the steady-state Boltzmann equation for Maxwell molecules appearing 
m eq. (gUJ), while 4>^s are different from each other. 



IV. THERMODYNAMIC QUANTITIES 

We can introduce the general form of the heat flux as 

J x = -zuT^d x T } (41) 

where (p indicates temperature dependence of the thermal conductivity and vj is a constant 
that depends upon microscopic models. J x is constant from the solubility condition (JBIJ). 
For example, tp is calculated as 1/2 for hard-core molecules and 1 for Maxwell molecules; w is 
determined as 756h/«(k / 'rnn) 1 1 2 / '64a 2 with bu — 1.025 for hard-core molecules ( see eq. (|35jl ) 
and 5k 2 (2m/ G) 1 / 2 / 4A for Maxwell molecules [lj], where G is the constant of proportionality 
between the intermolecular force and the reciprocal fifth power of the distance, and A is 



a number constant erroneously evaluated as 1.3682 by Maxwell 



12. 131 and recalculated as 



1.3700 by Chapman 43]. Note that ip and w cannot be determined explicitly from the BGK 
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equation. From eq. (j4~Tj) . the temperature profile T(x) in the nonequilibrium steady state 
can be determined as 

T(x) = [T(0)^ +1 -(^ + l)^x]^r 

zu 

- T(0)[1 " wT(oy+i x ~ ^f^" 2 '' (42) 

to second order. The temperature profile T(x) becomes nonlinear except for ip = 0. This 
fact is in accord with a numerical result that a temperature profile becomes nonlinear as the 
heat flux becomes larger, j^J 

Using eq.(j2H), the pressure tensor Py in the nonequilibrium steady state can be obtained 

as 

to 

Pij = nnT[5 l] + ^-^i (43) 

with the unit tensor 5^ and the tensor components Ap given in Table IIV1 Note that the 
off-diagonal components of Ap are zero and X V p = Ap is satisfied. The values of Ap for 7th 
approximation b\ r , bo r and b 2r for hard-core molecules, i.e. 7th approximation Ap seems to 
be converged to three significant figures, as can be seen from Table IIV1 We find that Ap 



for hard-core molecules differs from that for Maxwell molecules not only quantitatively but 
also qualitatively: Af? ^ X V p = Ap for hard-core molecules in the nonequilibrium steady 
state, while \p x = \ V p = Ap 2 = for Maxwell molecules. For hard-core molecules, P xx 
becomes smaller than P yy and P zz regardless of the sign of J x . It is also important that, 
since nnT is required to be uniform from the solubility conditions for <ffi*' in eq.(|S]l. we find 
from eq. (|42jl that in eq. (|4Hjl is uniform to second order in the steady state. Additionally, 
since the temperature profile T(x) has already been given in eq.()42jl. we can determine the 
density profile n(x) to second order for uniform P^. Note that the equation of state in the 
nonequilibrium steady state is not modified to first order, and that Ap for the steady-state 
BGK equation is identical with that for Maxwell molecules. 

Each component of the kinetic temperature in the nonequilibrium steady state, i.e. 7] 
for i = x, y and z is calculated as 

dc^— / 



771 / / 2 



ukT mJl , , AA . 

— [l + ^T*M> (44) 
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for i = x, y and z. Values for the constants At; for i = x, y and z are given in Table ITVl Note 
the relation At„ = Xt z - At; with 7th approximation b\ r , b 0r and bi r for hard-core molecules, 
i.e. 7th approximation At; seems to converge to three significant figures, according to the 
results reported in Table II VI Since tlkT and J x are uniform, we see that the correction 
term of T x in eq . (|44|) is independent of x. The value of At- for hard-core molecules is 
seen to be qualitatively different from that for Maxwell molecules: \t x 7^ Xt v = At z for 
hard-core molecules in the nonequilibrium steady-state, while At, = \t v = At z = for 
Maxwell molecules in that state. For hard-core molecules, T x becomes smaller than T y and 
T z regardless of the sign of J x , which means that the motion of hard-core molecules along the 
heat flux becomes dull. We note that 7] for the steady-state BGK equation is identical with 
that for Maxwell molecules, and that, to first order, 7] for hard-core molecules is isotropic. 
The Shannon entropy in the nonequilibrium steady state S is defined via 



Values for the constant As are given in Table IIVI As for 7th approximation b\ r , b 0r and 
6 2r - for hard-core molecules, i.e. 7th approximation As seems to converge to four significant 
figures, see Table UVl It is found that As for hard-core molecules is close to that for Maxwell 
molecules, while the latter is identical with that for the steady-state BGK equation. Note 
that the Shannon entropy in the nonequilibrium steady state is not modified to first order. 

V. DISCUSSION 

It is noteworthy that all macroscopic quantities now become able to be calculated directly 
from the explicit velocity distribution functions of the steady-state Boltzmann equation to 
second order. Actually, as in Appendix |Ej we examine the nonequilibrium steady-state ther- 
modynamics(SST) proposed by Oono and Paniconi|45(], and extended by Sasa and Tasaki 
by using the velocity distribution function to second order for hard-core molecules we have 
derived in §1111 and that for Maxwell molecules. As a result, our results do not support 
SST. We mention that the full information of the second-order coefficients B^ r and B l ^ r is 
necessary for the calculation of some physical quantities in the test of SST. (see Appendix 
IE]) Furthermore, as indicated in the introduction, the other important application of the 



3 




(45) 
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explicit second-order velocity distribution function for hard-core molecules is to evaluate the 
nonequilibrium effect on rate of chemical reaction of gases. The line-of-centers model j^] 
which is accepted as a standard one requires the full form of the second-order velocity distri- 
bution function for hard-core molecules j^J. Our result enables us to calculate the reaction 
rate.U 

We have found that there are qualitative differences between hard-core molecules and 
Maxwell molecules in the nonequilibrium steady state: second-order corrections appear for 
hard-core molecules in the pressure tensor Pij and the kinetic temperature 7], while no 
correction to these quantities appears for Maxwell molecules, as Table HY1 shows. It should 
be noted that the qualitative differences between hard-core molecules and Maxwell molecules 
still appear no matter which boundary condition is adopted, that is, the isotropy and the 
anisotropy of the pressure tensor in eq. (|43jl and the kinetic temperature in eq.([44p are not 
affected by any kinds of boundary conditions. It is conjectured that these differences are 
attributed to the special nature of Maxwell molecules: Ak r for Maxwell molecules, i.e. Aj^. 
is independent of the magnitude of the relative velocity g, because F£ r (x) is independent of 
g. In general, Ak r for molecules which interact with each other by a central force depends 
on g, because F£ r (x) generally depends on g. Therefore, it may be suggested that molecules 
which interact via a central potential still have the qualitative differences from Maxwell 
molecules: second-order corrections in pressure tensor and kinetic temperature 7] may 
appear also for such molecules as well as hard-core molecules. 

It is also found that the pressure tensor P^ and the kinetic temperature 7j for the steady- 
state BGK equation are qualitatively different from those for the steady-state Boltzmann 
equation for hard-core molecules but agree with those for the steady-state Boltzmann equa- 
tion for Maxwell molecules, as illustrated in Table ITVl Since the Chapman- Enskog solution 

of the steady-state BGK equation is asymptotically correct j^. y, y, y, q 



we may 



conclude that the steady-state BGK equation does not capture the essence of hard-core 
molecules, but captures that of Maxwell-type molecules. This conclusion indicates the pos- 
sibility that even the exact solution of the steady-state BGK equation can be applied only 
to Maxwell-type molecules. 

Finally, we consider the possibility of the existence of a universal velocity distribution 
function in the nonequilibrium steady state. When we express the velocity distribution 
function for hard-core molecules to second order using the heat flux J x as in eq. (|39jl . it be- 
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comes independent of the diameter a of the hard-core molecules. However, the explicit form 
of the velocity distribution function for hard-core molecules f!39p definitely differs from the 
precise form of the velocity distribution function for Maxwell molecules (|4U|) as FigEI shows. 
Actually, as mentioned above, we have shown that the results calculated from the former 
are qualitatively different from those calculated from the latter. These results indicate that 
the characteristics of microscopic models appear in a nonequilibrium steady state solution 
and affect even qualitatively the macroscopic quantities of a system in that nonequilibrium 
steady state. It can be concluded that a universal velocity distribution function does not 
exist in the nonequilibrium steady state even if the velocity distribution function is expressed 
only in terms of macroscopic quantities. [5] 



VI. CONCLUSION 



The velocity distribution function of the steady-state Boltzmann equation for hard-core 
molecules subject to a temperature gradient has been derived explicitly to second order 
in the density and the temperature gradient as shown explicitly in eq.()39p and Fig|21 In 
the nonequilibrium steady-state system, qualitative differences between hard-core molecules 
and Maxwell molecules are found in the pressure tensor ()43j) and the kinetic temperature 
f)44jl : it appears that Maxwell molecules do not possess the characteristics of other models 
of molecules which interact with each other by central forces in the nonequilibrium steady- 
state system. Additionally, we have found that the steady-state BGK equation belongs to 
the same universality class as Maxwell molecules, and that it does not capture the essence of 
hard-core molecules. We finally conclude that a universal velocity distribution function does 
not seem to exist, as Fig|3] explicitly shows, even when the velocity distribution function is 
expressed only in terms of macroscopic quantities. 
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APPENDIX A: THE CALCULATION OF Vt kr . 

From the definition of Qkri ^kr can be calculated using t 
the spherical harmonic functions and the Sonine polynomials. 
Ek r and Gk r defined as 



re mathematical properties of 



One can calculate Dk r , 



1 1 1 

Dkr — *^ ClxQkr ^* avt Efo r = <C ClyQkr ^av and Gkr — ^ ClzQk r ^> a v • (Al) 

n n n 

The results can be written as 

D k,r = + r + l)B^ 1>r - S« 1|P _J - ^j[B^,r ~ (r + l)B { k \ r+1 ], (A2) 

and 



2fc + 3 1 



and 



Gk ' r = WTZ^ + T + l^ Bk+1 < r ~ B k+hr-l] - 2k^l^ Bk ~ llT " ^ + l ) B *-l.r+l]- ( A4 ) 

One can also obtain 

d T k 

< c lx d x Q kr >av= -n-7p-[(r + ^)D k , r - D k , r -i}- (A5) 

Similarly < c\ y d y Qk r >av and < c\ z d z Qk r >av are obtained by replacing the differential 
coefficients with respect to x by the corresponding differential coefficients with respect to 
y and z, the D^s by the corresponding E^s and Gfc, r 's, respectively. Substituting these 
results into eq.(j2IJ), &kr finally becomes eq.(J2EJ)- 
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APPENDIX B: THE DETAILS OF F^fo) 



The details of F^ r (x) are as follows. Substituting the general forms of /i, f 2 in eq.flHJ) 
and Q' fcr in eq.(|22J) into F^ r (x) in eq. (JZ3j) . F^.(x) can be written asfl^ 

fUx) = e ^r^s^i 2 '' 1 '' 2 ^), (bi) 

ni,ri2,fci,fc2,fci>Pi>0,fc2>P2>0 

where the summation with respect to pi and p 2 is performed from to k± and k 2 , respectively 
as seen in eq. ffTB|) . Here S^,'^'^' Pl ' P2 (x) is t ne characteristic integral defined as 



x 



exp[-(c? + ^)]^ +i (cf)5 fc " i 1 +i (c?)^ +i (c2)^dedc 2 dc 1 . (B2) 



2 n.i-r 5 «Tj 



Note that the integrals containing products like Yjz (ci)Z^ 2 (c 2 ) are zero, owing to 
the orthogonality properties of the spherical harmonic functions, while those containing 
Z^\ci)Z^\c 2 ) are identical with the corresponding integrals H£^ 1 '^' Pl ' P2 (x)- The factor 
W fe ™f ,P2 in eq.jEH) is defined as 

iT/Pi>P2/ d(pi) r>(pa) , /^f(pi) /~r(pa) \ / R q\ 

which is obtained from the prefactors and the coefficients in the general form of fi, f 2 
in eq.(JHJ) and Q' kr in eq.(j22J)- Note that appear from the integrals containing 

Z%°(ci)z£\c 2 ), and that B$ = B kr and C^ = from eqs.fllEJ), (HZj), JEJ and JED. The 

constant ^f 1 '^ is defined as 

*EJ = 1 for Pl =p 2 = 0, (B4) 

and 

Here ^f^ 2 for pi 7^ p 2 is not necessary for our calculation. Q It is found that we 
need only to evaluate the characteristic integral S^_^ 2,pi ' pi (x) i n order to calculate F^ r (x). 



Our calculation of 2^.'^'^' pi ' pi (x) is written in ref. [49(. Our calculation has been performed 
mainly based on Burnett's met hod 1-Y. We have, however, made some modifications on his 
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method, which make the calculation of H£^'^' Pl ' Pl (x) much easier. [49| We emphasize that 
our calculation could not be carried out completely if we did not make the modifications on 
Burnett's method. 

Once the characteristic integral S^.'^'^' Pl ' Pl (x) has been derived, F kr (x) is now calcu- 
lated from eq.flEU) with W^ 2 k f im in eo.l|B5|). Note that the coefficient term B^B^ + 
CUC^I in WJX ,fcf 1,P1 can be determined for each order when the suffices kx, k 2 , ni, n 2 , 
Px and p 2 are specified. 



APPENDIX C: THE CALCULATION OF THE FIRST ORDER COEFFICIENTS 

B kr 

Let us explain how to obtain the first-order coefficients, that is, how to solve the integral 
equation (JHJ). To begin with, we calculate in eq. (j26j) to first order; for first order 
corresponds to the right-hand side of eq.(jHJ). For first order, in eq. (|2fjjl can be calculated 
by substituting B 00 = 1 into the expressions of D k ri E k ^ r and Gk, r in eqs. f)A2|) . ()A3|) and 
flA4|) : the coefficient B 00 = 1 corresponds to fx = f[°\ and higher-order terms do not appear 
in for first order. for first order finally becomes 

lr T \2kTJ dz ' V ; 

Now f2j? r for first order is found to vanish unless k = 1, so that we need calculate only 
A^ r for first order; as was mentioned in the end of j^nj we do not need to consider the case in 
which the right-hand side of eq.© becomes zero.|2] To derive A^ r in eq. (|28jl for first order, 
we must calculate both W r 1 n ^ n ^ Pl ' Pl and 2^^'^' Pl ' Pl in F^ r (x) of eq. ljBfj) for first order as 
was shown in the previous subsection. The result for A^ r to first order can be written finally 
in the form 

^l = BooJ2 B l 1 K:if A °, (C2) 

"1 

where the set of the coefficients B\ ni B 00 is obtained from q' '° in eq. ()B3|) . fx in eq.pUjl 
contains only Bqq = 1 and the first-order coefficients, i.e. the family of B\ to first order; 
f 2 in eq. (j2T?j) also contains B 00 = 1 and the family of B\ to first order. Thus, we obtain 
only the term B\ ni B 00 from W^l' ' '° to first order using the fact that F kr (x) = unless 
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k = \ki — A^l + 2q. Note that it is sufficient to consider only the case for k\ > k2 as is 



explained in refs. 



HQ 



i,jr,ni, 0,0,0 
iw 1,1,0 



and that we set g = 0. The matrix M 1 ' 1 1 : 

3n 2 cr 2 m 4 rii! 
647tIk 4 T 4 



r,m, 0,0,0 



is thus obtained 



x 



r ,— 'T\ni, 0,0,0/ \ ^r,ni,0,0,0/nM X X j 

Fi,i,o (X) - "i,i,d (°)] sm 2 cos g d *' 



(C3) 



using eqs. 

For k — 1, eq. 
B{ ni , i.e. 



Bll) and (IB3I). 



gives simultaneous equations determining the first-order coefficients 



0,0,0 
J 



(C4) 



from eqs. ()Cl|) and ()C2|) . Note that we need only to obtain the first-order coefficients B\ ni 
for rii > 1, because B w = from eq. (|2(J|) . We have calculated the matrix M['[' 



r,m ,0,0,0 



.0 



for 



1 < r < 7 and 1 < ni < 7 from eq. ()C3)) . and we have confirmed that M^y ' 0,0 for 



1 < ^i < 7 calculated from eq. ()C3|) vanishes. Our result for M['"q ' ' for 1 < r < 7 and 
1 < ni < 7 is given in ref.j^]. At last, we can determine the first-order coefficients B\ by 
solving the simultaneous equations (|C4|) . that is, B\ ni can be obtained as 

r>l 



(C5) 



where X -1 represents the inverse matrix of a matrix X. Note that we have confirmed, using 
eqs. (jCl|) and (|C3|) . that both sides of eq. (|C4|) for r = vanish, so that we need only calculate 
both sides of eq. (jC4|) for r > 1. Finally, the results of the first-order coefficients B\ , i.e. 
B\ T in eq. (fTH|) can be calculated as in eq.(|HH|). 



APPENDIX D: THE CALCULATION OF THE SECOND ORDER COEFFI- 
CIENTS Bj} r 

Let us explain how to obtain the second-order coefficients, that is, how to solve the integral 
equation ©. The coefficients of first order, i.e. the family of B\ r , have been obtained as 
given in eq. (|33)) . so that we can employ them to determine the second-order coefficients, i.e. 
the family of B^ r . 

To begin with, we calculate fl^ r in eq. (j2Tj) for second order; for second order corre- 
sponds to the first term on the right-hand side of eq.©. For second order, fl^ r in eq.tj^HJ) 
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can be calculated by substituting the family of B\ r into the expressions of D^ r , Ek <r and 
Gk, r in eqs. (jA2|) . ()A3|) and ()A4|) : other terms do not appear in for second order. The 
results of the tedious calculation of to second order finally become as follows. For k = 0, 



becomes 



for r = and 1, 



for r = 2, and 



ng. = o, (di) 



(r 2 + ^ - -)6i r - (2r - -)&i, r _ a + &i, r _ 2 



(D3) 



for r > 3. Note that values for the constants &i r are summarized in Table |U For k = 2, fi^, 
becomes 



and 



and 



«2o = 0, (D4) 



n| = -| v^r» (^)~* [2 {dzT)2 " (4T)2 " (w2] (&12 ~ M ' (D5) 



1 ( m \~2 
16 V^kct 2 T 2 

x [-r(r + l)6 lir+1 + 2rb lr - b x ^ x \ , (D6) 
for r > 2. For = 1 and k > 3, we find that fi^ r becomes 

n£ = 0, (D7) 

for any value of r. 

Next let us calculate Aj^, in eq. (J25|) for second order. In order to derive Af r for second 
order, we have to calculate W^™lf 1,Pl and in F£ r (x) of eq. ()Bl|) to second order, 

as was shown in Appendix El For even k, Aj? r to second order results in the general form: 



aH _ r \^ R II » # r,ni,0,0,0 . \ ^ R I R I , # r,ni,n 2 ,0,0 

^kr - A)0 2^ ^k-2q, ni M k,k-2q,0 + J °lni J °ln2 M fc,l,l 

™i,0<<?<§ n i' n 2 

+ ^[B^BSI + C^C^WIX^ 1 ' 1 - (D8) 



ni,n 2 
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Here the set of the coefficients Bj}_ 2qni B 00 can be obtained from W^]^ 2q \ in eq. (jB3|) . The 
set of the coefficients B\ B\ is obtained from W^™ 2,0 ' m eq. (|B3J) : B\ from f\ and 
B\ n2 from f 2 are the first-order coefficients obtained in eq. (j3*3*Jl . so that B\ ni B\ n2 is second 
order. Similarly, the set of the terms B^B^ + C[^C[^ 2 is obtained from W^ 1 ^ 2 ' 1 ' 1 in 
eq. (|B3|) ; B^, from j\ and B^, C[^ 2 from f 2 are the first-order coefficients obtained 
in 1111 C 11 so that B^B^ and C[^C[]^ 2 are also second order. To second order, /i of 
eq.(j29j) contains only B Q0 = 1, the family of B\ obtained as eq. (jHH|) and the family of 

to be determined here; f 2 of eg. ((221) also contains only B 0Q = 1, the family of B\ 
and the family of B^ n to second order. Thus, we can only obtain the sets of the terms in 
eq. (|D8|) for second order by using the fact that F£ r (x) = unless k — \ki — k 2 \ + 2q: the 
second and the third terms on the right-hand side of ea. (ID8l) do not appear for odd k. Note 

nn 

that it is sufficient to consider only the case for ki > k 2l as is explained in refs.|15L I49|. and 
that B^c[^ 2 or c[^B^ 2 does not appear owing to the orthogonality properties of the 
spherical harmonic functions, (see Appendix |Bj) The matrix M^'^I^q in eq. (jD8|) is thus 
obtained as 



M, 



•r,ni,0,O,0 71 ° ( 771 V f 771 \ k ~ Q WU , ^ r 



x J\^°J(X) ~ S££S£f(0)] sin | cos |d X , (D9) 



using eqs.(|2HI), (|BT|l and (|E3|). Similarly, the matrices M^\' n2 ' u ' u and M^i in eq. 
are derived as 



x \\mr^\x) - *zt A0 (°)] sin f c ° s f d *> (dio) 



and 

2 \2tikTJ \2kTJ v 2 ;v 

x ^ V^i n2,1,1 (x) - H^r' 1 '^ )] sin f c ° s f d ^ ( Dn ) 

respectively. 

For even k, eq. (j32j) finally leads to simultaneous equations to determine the second-order 



23 



coefficients B^_ 2qni : 

nH _ \^ oil ,^r,ni,0,0,0 \^ R I R I , # r,ni,n 2 ,0,0 

S 'fer — n k-2q,n 1 1V1 k,k~-2q,0 ^ ^lni ^lna iW fc,l,l 

«i,0<g<-§ ™i- n 2 

+ E^n^S + ^^fflKS'" 2 ' 1 ' 1 , (D12) 
m,n2 

from eq. ()D8|) . The second and the third terms on the right-hand side of eq. (JD12[) correspond 
to J(fx , f^) in the integral equation (jOJ). Thus, it had been believed that eq. (|D12J) should 
be considered for all even k because the contribution from the right-hand side of eq.flHJ) would 
not become zero even when which corresponds to the first term on the right-hand side 
of eq.© is zero.Q] However, for k = 4, 6 and 8, we have confirmed that the second and 
the third terms on the right hand side of eq. (|D12|) disappear, which leads to the fact that it 
is not necessary to consider eq. (|D12|) for k = 4, 6 and 8. Therefore, it is natural to expect 
that we do not need to consider eq. (|D12J) for even k furthermore. Our contention that the 
second and the third terms on the right hand side of eq. (|D12|) will disappear for k = 2n 
with the integer n > 2 should be demonstrated by a mathematical proof in the future. It 
should be mentioned that the second and the third terms on the right-hand side of eq. (|D12j) 
do not appear for odd k, and to second order is found to be zero for odd k, so that 
we do not need to calculate Aj? r for odd k; it is not necessary to consider the case in which 
contribution from the right-hand side of eq.(jBJ) becomes zero.Q] 

Now, we need to consider eq. (|D12J) only for k = and 2. If k = 0, eq. (|D12)) leads to 
simultaneous equations to determine the second-order coefficients B^ ni : 



nH _ R II ji^nm ,0,0,0 , R I R I * # r,ni,n2,0,0 

"Or - 2^ if Oni iM 0,0,0 + 2^ B ln 1 B ln 2 M 0,l ) l 

ni>2 ni,n2 

+ E^^+^^SKrr- 1 ' 1 ' ( Di 3) 



711,712 



that is, 



R II _ J O h R 1 R 1 A/f"'™ 1 '™ 2 - ' 

B 0ni - Z^) 0r ~ ±, lni- fc »ln a M 0,l,l 

r>2 \ 7ii,?i2 
ni,n 2 J 

with Q^ r in eqs. (|Dl|) . (|D2|) and ()D3|) . We should derive the second-order coefficients B^ ni 
only for n\ > 2, because B 00 = 1 and B i = from eq.(|2~Uj). We have calculated the 
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matrix Mq'qq' ' ' for 2 < r < 6 and 2 < n\ < 6 from eq. (jD9j) . and we have confirmed 
that Mq'qq 0,0,0 vanishes for r = 0, 1 and 2 < n± < 6. We have also calculated the matrices 
MQ , ™ 1 1 ,ra2,0,0 for 2<r<6, l<n 1 <7 and 1 < n 2 < 7 from eq. ()D10|) . and we have 
confirmed Mq'™ 1 ^" 2,0,0 vanishes for r = 0, 1, 1 < n\ < 7 and 1 < n 2 < 7. Our results for 



Koo ' ' for 2 < r < 6 and 2 < ni < 6 and M^ n2 ' u ' u for 2 < r < 6, 1 < n x < 7 and 



r,m,n2,0,0 



1 < ^2 < 7 are given in ref. 42|. The matrices Mg'™ 1 / 12 ' 1 ' 1 for 2 < r < 6, 1 < rii < 7 and 
1 < "^2 < 7, which can be calculated from eq. (|Dll|) . are also obtained from the symmetric 
relation M^™' ' = M r , " 1 1 ' n2 ' ' arising from properties of the spherical harmonic function. 
Finally, we can determine the second-order coefficients B^ in fi, i.e. B 1 q t in eq.(JTlJ) as in 
eq-P7|l. 

If k — 2, eq. ()D12|) leads to simultaneous equations to determine the second-order coeffi- 
cients i?2ra : 



ni>0 ni,n2 

+ E^^+^^ffl^xr' 1 ' 1 ' ( Dl5 ) 



ni,n 2 



that is, 



< = E K - E ^n^L.^r- ' 

r>0 v. ni,ri2 



^ [B^X + C^lC^M^ 2 ' 1 ' 1 \ (M^ 0,0,0 )- 1 , (D16) 

m,n2 



with Qf r in eqs. (jD4j) . (fD5j) . and (jDEJ). Note that we have confirmed M^' 0,0 ' in eq. (pTl5jl 
becomes zero, which had been also confirmed by Burnett 15]. We should derive the second- 
order coefficients B% for rii > 0. We have calculated the matrix ' ' ' for < r < 6 
and < ri\ < 6 from eq. (jD9|) . and also the matrices M^'^'" -2 ' ' for 0<r<6,l<ni<7 
and 1 < ri2 < 7 from eq. (|D10|) . Our results are given in ref. j^. The matrices M^'iV™ 2 ' 1,1 
for < r < 6, 1 <■«,!< 7 and 1 < n 2 < 7, which can be calculated from eq. (jDll|) . are also 
obtained from the symmetric relation Mj^'"' 2 ' ' = — 2M2'™ V" -2 ' 1 ' 1 arising from properties of 
the spherical harmonic function. [3] The second-order coefficients Bj} ini in fx, i.e. B l ^ r in 
eq.f|l(j|) can be written in the final form shown in eq.(|38|). 
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APPENDIX E: TEST OF THE NONEQUILIBRIUM STEADY-STATE THER- 
MODYNAMICS 



We examine SST from the microscopic viewpoint, that is, by applying the velocity distri- 
bution function for the steady-state Boltzmann equation for both hard-core molecules and 
Maxwell molecules to the cell on the right-hand side in a simple nonequilibrium steady-state 
system illustrated in FigJU We also use the velocity distribution function for the steady-state 
BGK equation to second order in the test of SST. This simple nonequilibrium steady-state 
system is inspired by SST suggested by Sasa and Tasaki[46]. They predicted the following 
results. When the system shown in Fig0] is in a steady-state, the osmosis, defined as the 
difference between the pressure P xx of the cell in the nonequilibrium steady state and the 
value Pq of the cell at equilibrium, namely 

AP EE P xx - P , (El) 

always becomes positive. Additionally, there is a relation 

n(0) (dP x 



n \dP J ToJx 

connecting n(0), P xx , no and Pq, where n(0) is the density of the cell in the nonequilibrium 
steady state around the hole and uq is that of the cell at equilibrium. 

We consider the nonequilibrium steady-state where the mean mass flux at the hole is 
zero: 



poo poo pO poo 

/ dc x / d 2 c±mc x f + / dc x / d 2 c±mc x f\ x= o = 0, 

J J— oo J —oo J —oo 



(E3) 



where c± represents the components of the velocity which are orthogonal to c x , i.e. c y 
and c z . Note that we consider the mean mass flux at the hole, so that we put x = in 
/■ fo = n (m/27rKT ) 3 / 2 e~ rnv2 / 2KTo is the velocity distribution function of the cell on the 
left-hand side at equilibrium at temperature T and density no- From eq. (jE3|) . the relation 
between n(0) and n is obtained as 

TflJ 2 

n(0)=n o [l+A n - 5 - 5 ^], (E4) 

to second order. The value for the constant A n is given in Table |3 The density of the cell 
in the nonequilibrium steady state around the hole n(0) is greater than that of the cell at 
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equilibrium no regardless of the sign of J x . We emphasize that A n could not be calculated if 
we did not derive the explicit form of the velocity distribution function / to second order. 
We have adopted the boundary condition around the hole, T(0) = T , in order to examine 
SST. Although it is believed that the Knudsen layer effect, i.e. the slip effect is dominant 
around the 'wall'0, 16, £|, for reasons of simplification, we do not consider the slip effect 
around the 'hole' in this paper. If the slip effect is not dominant, the density of the cell in 
the nonequilibrium steady state will always be larger than that of the cell at equilibrium, 
regardless of the sign of J x . This can be tested by experiments on the nonequilibrium 
steady- state system shown in FigJH 

We also calculate the osmosis AP in eq. (jEl|) as 



AP = (E5) 



to second order, using eqs.([43j). (|K4|) and P = n nT . Values for Aap are given in TableJ 
We have found that Aap is always positive, which agrees with the prediction by SST 46]. 
Furthermore, we are able to test the relation (jE2|) . Though substitution of eqs.([43|) and 
(|E4|) into eq. (|E2|) leads to the relation X x p /X n = —2, our numerical results conflict with this 
relation: our results give Af?/A n = —0.5604 for 7th approximation bo r and &2r for hard-core 
molecules and Af?/ A n = for both Maxwell molecules and the steady-state BGK equation. 
We have also confirmed that the relation Af?/A n = —2 predicted in is not modified 

if the boundary condition around the hole can be written in the form: 

to 

T(0)=T [1 + A T -^], (E6) 

where the constant Xt represents the difference between the temperature of molecules around 
the hole and that of the mid-wall. Our boundary condition T(0) = To corresponds to putting 
Xt = for any kinetic models. It can be concluded that, although we regard a state in which 
the mean mass flux at the hole is zero to be a nonequilibrium steady state, as in eq. (jE3|) . 
this state has yet to be interpreted phenomenologically. 

On the other hand, by virtue of the derivation of the explicit form of the velocity distri- 
bution function / to second order, we can also calculate the x component of the heat flux 
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at the hole as 

[2kT\ z f 00 , f°° . mc 2 r [2kT\ z f° f 00 l2 mc 2 „, 
J = / dc x / d ci— cj + / dc x / d c±— — c x f \ x=0 

to second order using eq. ()E4jl . Numerical values for the constant Aj* are given in Table 
fVl The first-order term on the right-hand side of eq. (jE7|) is from Fourier's law, while the 
second-order term also appears on the right-hand side of eq.([E7[). though the second-order 
heat flux j( 2 ) does not exist in the cell in the nonequilibrium steady state. This fact indicates 
that T(0) is not appropriate for the nonequilibrium temperature, if J* suggests the existence 
of the nonequilibrium temperature. 
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TABLE I: Numerical constants b\ r in ea. (|33j) 



r r < 4 r < 5 r < 6 r < 7 Maxwell's &i r 



1 1.025 


1.025 




1.025 




1.025 




1 


2 4.881 x 1(T 2 


4.889 x 10~ 


-2 


4.891 x 10- 


-2 


4.892 x 10- 


-2 





3 3.639 x 1(T 3 


3.698 x 10" 


-3 


3.711 x 10- 


-3 


3.715 x 10- 


-3 





4 2.526 x 10~ 4 


2.838 x 10- 


-4 


2.905 x 10- 


-4 


2.922 x 10- 


-4 





5 


1.855 x 10- 


-5 


2.123 x 10" 


-5 


2.187 x 10- 


-5 





6 






1.284 x 10- 


-6 


1.492 x 10- 


-6 





7 










8.322 x 10- 


-8 






TABLE II: Numerical constants &or in ea, ()37|) 



r 


r < 4 


r < 5 


r < 6 


r < 7 


Maxwell's &or 


2 


4.434 x 10" 1 


4.390 x lO" 1 


4.381 x 10" 1 


4.380 x 10" 1 


825 
1024 


3 


-4.935 x 10~ 2 


-5.342 x 10~ 2 


-5.413 x 10~ 2 


-5.429 x 10" 2 


25 
256 


4 




-3.581 x 10~ 3 


-4.007 x 10~ 3 


-4.098 x 10" 3 





5 






-2.779 x 10~ 4 


-3.184 x 10~ 4 





6 








-2.087 x 10~ 5 
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TABLE III: Numerical constants bi r in ea. (|38j) 



r 


r < 4 


r < 5 


r < 6 


r < 7 


Maxwell's &2r 





-3.353 x 1(T 2 


-3.327 x 10~ 2 


-3.322 x 10~ 2 


-3.320 x 10~ 2 





1 


-1.285 x KT 1 


-1.278 x 1CT 1 


-1.277 x 10 _1 


-1.276 x 10 _1 


75 
896 


2 


6.320 x 1CT 1 


6.394 x 10~ 2 


6.410 x 10~ 2 


6.414 x 10~ 2 


125 
1536 


3 


4.884 x 10~ 3 


5.395 x 1(T 3 


5.496 x 10~ 3 


5.521 x 10~ 3 





4 




3.609 x 10- 4 


4.101 x 10~ 4 


4.214 x 10" 4 





5 






2.685 x 10~ 5 


3.106 x 10~ 5 





6 








1.861 x 10~ 6 






TABLE IV: Numerical constants for the macroscopic quantities I: the ith approximation quantities 
for hard-core molecules and the exact values for Maxwell molecules and the steady-state BGK 
equation. 



ith 


Xf 




xf 








X T y 




Xs 




4th 


-4.647 x 10" 


-2 


2.324 x 10" 


-2 


-2.324 x 10" 


-2 


1.162 x 10" 


-2 


-2.034 x 10" 


-1 


5th 


-4.610 x 10" 


-2 


2.305 x 10" 


-2 


-2.305 x 10" 


-2 


1.153 x 10" 


-2 


-2.035 x 10" 


-1 


6th 


-4.602 x 10" 


-2 


2.301 x 10" 


-2 


-2.301 x 10" 


-2 


1.151 x 10" 


-2 


-2.035 x 10" 


-1 


7th 


-4.600 x 10" 


-2 


2.300 x 10" 


-2 


-2.300 x 10" 


-2 


1.150 x 10" 


-2 


-2.035 x 10" 


-1 


Maxwell 






















l 

5 




BGK equation 






















1 

5 
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TABLE V: Numerical constants for the macroscopic quantities II: the ith approximation quantities 
for hard-core molecules and the exact values for Maxwell molecules and the steady-state BGK 
equation. The values of A n , Aap and Aj* for hard-core molecules has not yet converged to 4th 
approximation &o and b<i r values. The ratios of the 7th to the 6th approximation A n , Aap and Aj* 
are 1.011, 1.023 and 1.004, respectively, so that the errors included in the 7th approximation A n , 
Aap and Aj* appear to be less than about one or two percent. Note that ra(0) = no to first order, 
and that osmosis does not appear to first order. 



ith 


An 






Aap 


Aj* 






4th 


9.255 x 


10" 


-2 


4.608 x 10~ 2 


-3.237 x 


10" 


-1 


5th 


8.528 x 


10" 


-2 


3.917 x 10~ 2 


-3.109 x 


10" 


-1 


6th 


8.296 x 


io- 


-2 


3.694 x 10~ 2 


-3.073 x 


io- 


-1 


7th 


8.210 x 


io- 


-2 


3.609 x 10~ 2 


-3.060 x 


10" 


-1 


Maxwell 


71 

157E 






71 
1575 


41 
105 






BGK equation 


2 

25 






2 

25 


11 

25 








FIG. 1: Schematic description of an interaction. 
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Cx 



FIG. 3: Comparison of the scaled (f>( 2 ^s for hard-core molecules with the scaled (j)^ for Maxwell 
molecules. The dashed line, the dash-dotted line, the long-dashed line and the solid line correspond 
to the scaled <f>^s for hard-core molecules with 4th, 5th, 6th and 7th approximation bo r s and &2rS, 
respectively. The dotted line is the scaled for Maxwell molecules. Note that we put c y = c z = 0. 
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FIG. 4: The cell on the left-hand side is at equilibrium at temperature To. The cell on the right- 
hand side is in a nonequilibrium state under a temperature gradient along the x-axis caused by 
the right wall (at x = Lq) at temperature T\ and the thin mid-wall (at x = 0) at temperature 
To, of thickness less than or equal to the mean free path I of the dilute gases. Both cells are filled 
with dilute gases and connected by a small hole of diameter d on the thin mid-wall. The diameter 
of the small hole d is much smaller than /, i.e. d <^ I. Molecules which have passed through the 
small hole are relaxed into the state of the cell they go into after a few interactions, so that they 
do not affect the macroscopic state of that cell. The mean mass flux at the small hole is zero in 
the nonequilibrium steady state. 



36 



